FINITE SPEED OF PROPAGATION AND LOCAL 
BOUNDARY CONDITIONS FOR WAVE EQUATIONS 
WITH POINT INTERACTIONS 
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Abstract. We show that the boundary conditions entering in the 
definition of the self-adjoint operator A A ' B describing the Lapla- 
ry | ■ cian plus a finite number of point interactions are local if and only 

if the corresponding wave equation <f> — A A,B (j> has finite speed of 
. propagation. 

t-H ! 

> ; 

CS| ' 1- Introduction and framework 

O '■ 

CN ■ The Laplace operator with point interactions in L 2 (M 3 ) can be de- 

fined in the following way (see p[J [2] and references therein for more 
q \ details). Consider any finite set Y = {%}™ =1 of points from M. 3 . Then 

the restriction — A of the Laplace operator —A to the set of functions 
Q_c from the Sobolev space H 2,2 (R 3 ) vanishing at all points yj 6 Y is a 

symmetric operator with deficiency indices (n,n). The domain of the 
adjoint operator — Aq coincides with H 2,2 (M. 3 \ Y). Every function <f> 
from H 2 ' 2 (M. 3 \ Y) possesses the following asymptotic representation in 
a neighborhood of each point yj E Y 

><: <K*) = T- n -^ s J +¥ J + 0{\x-y ] \), x^y v (1) 

Att\x — yj\ ■> ■> 

where the coefficients (jfj and $ can be considered as certain general- 
ized boundary values of the function (p. These boundary values can be 
used to describe all self-adjoint extensions of — A as restrictions of the 
adjoint operator — Aq to the set of functions satisfying the generalized 
boundary conditions: 

Af = Bcf s . (2) 

The n-dimensional complex vectors <p s and <p r have coordinates <pj and 
<f) r j respectively and the n x n matrices A, B satisfy 

AB* = BA\ (3) 

and 

rank (A, B)=n, (4) 
l 
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where (A, B) is the nx2n matrix obtained from A and B. The boundary 
form for the operator — Aq 



((-A$)0,V)l» - fa, (-AJ)V> 



L- 



determines the following symplectic form in the space C 2n 9 = 
(0 s , 0' r ) of boundary values: 



0,-0 



I n \( r 



(5) 

C 2n 



Then all self-adjoint extensions of — Ao can be described by Lagrangian 
planes associated with this symplectic form. Every such plane C is 
described by (J2J provided that the matrices A, B satisfy conditions 
(J3J) and (jlj). The first condition guarantees that the symplectic form 
vanishes for all 0, ip G C, i.e. that the corresponding extension of — Ao 
is symmetric. The second condition guarantees that the plane L has 
the maximal dimension n, i.e. that the operator extension is not only 
symmetric but self-adjoint as well. 

Definition 1. The operator — A A,B is the restriction of the adjoint 
operator — Aq to the domain 

D (-A A ' B ) = H 2 ' 2 (R 3 \ Y) : A$ r = B$ s j . (6) 

It is clear that different pairs of matrices (A, B) can determine the 
same self-adjoint extension of — Ao- If the matrix A is invertible then 
the boundary conditions (0) can be written as 

¥ = Hft, (7) 

with a Hermitian matrix H = A~ l B. The operator — A^ will be iden- 
tified with the operator — A In,H . 

The action of the operator — A A ' B coincides with the action of the 
(differential) Laplace operator and therefore this operator is always 
local, i.e. supp — A A)B (f) C supp0 for any <fi G D (—A A ' B ) and for any 
admissible pair A, B. However the boundary conditions (J2J) entering in 
Definition [T] are local, i.e. do not connect the "boundary values" 4>j,4> r j 
at different points x/j if and only if the matrices A and B can be chosen 
diagonal. In this case to check whether the function satisfies the 
boundary conditions (JJJ) or not, one needs to check the behavior of the 
function in a certain small neighborhood of each point yj separately. 
In this case the n-dimensional Lagrangian plane £ corresponding to (J2J) 
can be presented as a tensor product of n 1-dimensional Lagrangian 
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planes for the symplectic form in C 2 



o i \ i Vi 



-1 ) \ ib s , . 



defined for the boundary values (0£, 0|) associated with the point j/j. 
Such boundary conditions and corresponding self-adjoint extensions of 
— A will be called local. 

In fact all interesting local extensions corresponds to boundary con- 
ditions that can be written in the form ((7|) with a diagonal matrix H. 
The boundary conditions are local if the matrices A, B can be chosen 
diagonal. These conditions can be written as ((Zj) if the matrix A is 
invertible which can be assumed without loss of generality. Indeed if it 
is not the case then one can remove few points from the set Y to get 
boundary conditions equivalent to © with invertible matrix A (may 
be having a certain smaller dimension). 

The purpose of this note is to show (see Theorem 0] below) that the 
boundary conditions © are local if and only if the wave equation 

- A A > B = (9) 

has finite speed of propagation. We remind that in general an abstract 
wave equations 

+ „40 = O (10) 

is said to have finite speed of propagation if for any solution and for 
any real t, 

diam(supp <p(t j) < 2v \t\ + diam ^supp 0(0) U supp 0(0) j (11) 

holds with a certain v £ (0, oo). While it is obvious what "local bound- 
ary conditions" means in the case of — A A,B , the situation is not so clear 
for self-adjoint operators arising from extensions of symmetric opera- 
tors obtained by restricting the Laplacian to the set of smooth functions 
with compact support on R 3 \ M, where M is not discrete, for example 
a subset with Hausdorff dimension < d < 2 or a low dimensional 
manifold [21 IE] • Since the wave equation (fTT)|) is always well defined (in 
the sense that it generates a strongly continuous group of evolution) for 
any bounded from below self-adjoint operator A (see e.g chapter 
2, section 7), our result could be used to shed light on the problem of 
the locality of boundary conditions in more complicated situations. 

We conclude by pointing out that locality of boundary conditions, 
hence finite speed of propagation, is equivalent to locality in the sense 
of forms. According to [Hj, section 1.2, the self-adjoint operator — A A,B 
is said to be form-local if the sesquilinear form of the operator —A A,B 
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vanishes on any two functions from the domain D(F A ' B ) of the qua- 
dratic form having disjoint supports. The Theorem 2 from [7] implies 
that the operator — A ,B is form local if and only if the boundary 
conditions (0) can be written in the form (J7J) with local (i.e. diagonal) 
operator H. Hence the boundary conditions (J7J) with diagonal matrices 
H describe all form-local operators. 

2. Proofs 

The following representation for functions from the domain of the 
operator — A A ' B will be used in our proofs. 

We denote by H 1,2 (R 3 ) the homogeneous Sobolev space of tempered 
distributions such that iV0 is square integrable. Then H 2,2 (R 3 ) de- 
notes the space of G f/~ 1,2 (IR 3 ) such that — A0 G L 2 (M 3 ). In general 
neither G H 1 ' 2 ^ 3 ) nor e_H 2 ' 2 (R 3 ) imply G L 2 _(IR 3 ); by Sobolev 
embedding theorems one has i? 1,2 (R 3 ) C L 6 (M 3 ) andi? 2 ' 2 (R 3 ) C C b (R 3 ) 

We use the symbol Q to denote the Green function of —A, i.e. 

six) :=^n- 

An \x\ 

For any yj G Y we define the functions 

Gj(x) := Q{x - yj) , dj(x) := \x - yj\ , 
and the symmetric matrices G = (Gij) and D = (dij), 1 < i, j < n, by 
G i:j := Q{yi - yj) , i^j, Gjj := , d {j := \y { - yj\ . 

Lemma 2. The self-adjoint operator — A )B given in Definition 1 can 
be re-written in the following way 

D(—A A,B ) := J G L 2 (R 3 ) : = 0o + £ $ 0o G # 2 ' 2 (M 3 ), 

I l<j<n 



6 g c«, ^ A^j) = j2 ( B - AG h a I 

l<j<n l<J<n J 



-A^ B := -A0 O . 

Proof. It is sufficient to note that for any from the operator domain 
defined above one has G if 2 ' 2 (R 3 \ Y) and 

Jim (0-C^,)(x) = o ( % )+ £ G^Cj. 

Thus 



FINITE SPEED OF PROPAGATION AND LOCAL BOUNDARY CONDITIONS 5 



□ 

The proof of our main result relies on the following representation 
of the solutions of the Cauchy problem given by the wave equation 
(2). This result, in the simpler case of a single point interaction, was 
already obtained, by a different, less explicit proof, in (also see [§]). 

Theorem 3. The Cauchy problem 

= A A ' B <j> 
0(0) G D(A A ' B ) (12) 

0(0) g d\f a ' ) 

has an unique strong solution 

G C(R; D(A A ' B ) n C\R; D(F A ' B )) n C 2 (R; L 2 (R 3 )) 
explicitly given, when t > 0, by 

0(t) = 0/(t) + Yl W-di)<2ft-di)Gi > 

l<i<n 

where 9 denotes the Heaviside function, <fif is the unique solution of the 
Cauchy problem for the free wave equation 

h = A h (13) 

0/(0) = 0(0) (14) 
0/(0) = 0(0), (15) 

and C<j>{t)> t > 0, is the unique solution of the Cauchy problem for the 
system of inhomogeneous retarded first-order differential equations 

- oAuj) J + B^i = Yl A ^k0(- - d jk ) C{(. - d jk ) , 

/ l<J,fc<n 

(16) 

O(0) = C m ■ (17) 

Proof. Being — A A ' B a bounded from below self-adjoint operator, it is 
well known from the theory of abstract wave equations in Hilbert spaces 
(see e.g chapter 2, section 7) that the corresponding Cauchy prob- 
lem has an unique strong solution in C(R; D(A A ' B )nC 1 (R; D(F A ' B )) n 
C 2 (R; L 2 (R 3 )). Let us denote by <p(t) such a solution. By the structure 
of D{— A A ' B ) given in Lemma 121 we know that 

0(t) = o (t)+ Yl 

l<i<™ 
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with cj) (t) e H 2 ' 2 {M. 3 ) and 

E Aj Mt, vj) = E ( s " AG ^ ( 18 ) 

for all tel. Let us now define, for any t > 0, 

l<i<« 

where </>j(i) denotes the spherical wave 

^•(t) := eit-d^Qit-d^g,. 



Since 
one has 



0/ =a (0o - E & ) - E c}^ = a (00 - E - 

E <M = A0/ 

and thus 0/ is the unique solution of the Cauchy problem (fT3J) . (fTljl . 
(fTHJ) . Writing the boundary conditions (fTHjl by using the decomposi- 
tions of cf) given by the relation one obtains 

E - A ifc G fcj ) c}(*) 

l<7,fc<n 

= E 4* (ji m E (*) ) 

l<i<" y l<fc<n / / 

= E 4* ( E 9(t-d jk )G jk (i(t-d jk ) 

l<i<n V l<fc<n 

c^-i%--^i)-aw' 



E G iJfc Cj(t) + lim 



i< fc <„ ^ 47r l% - x \ 

<j<n \ l<fc<n 



1 ,, 



Kfc<n 
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This shows that ^ is a solution of (|TB|). (fTTj) . Such a Cauchy problem 
has an unique solution. Indeed let £(i), with £(0) = 0, solve the system 

By dHJ) and Ker( J 4)flKer(i?) = {0}, which is a consequence of (jlj), one 

has 

I det(iA + B)\ 2 = det((iA + B)(-iA* + B*)) = det(AA* + BB*) ^ 0, 

so that zA + B, z G C, is a regular pencil of matrices, i.e. : = 

det(zA + 5) is not the zero polynomial. By Theorem 3.2 of [3] this 
implies ((t) =0. □ 

The above theorem has an analogous version for negative times. In 
this case one obtains a representations involving, instead of the retarded 
waves 6{t — dj) — dj) Qj, the advanced ones 9(— t — dj) Q(i + dj) Qj. 
Since there is no substantial difference between these two situations, 
in the proof of the following theorem we will consider only the positive 
time case, the proof in the negative time case being essentially the 
same. 

The previous theorem shows that compactly supported Cauchy data 
always give rise to compactly supported solution. However this does 
not necessarily implies a finite speed of propagation. Indeed we have 
the following 

Theorem 4. Let —A A,B be a point perturbation of the Laplace operator 
described by Definition^ with the matrices A, B satisfying and Q). 
Then the wave equation 

4> - A A > B <f) = (20) 

has finite velocity of propagation if and only if — A A ' B is defined by local 
boundary conditions, i.e. if and only if both matrices A and B can be 
chosen diagonal simultaneously. 

Proof. Let (j>{t) be a solution of (fT2"j) with diam(So) < +00, where 
5*o := supp (f)(0) U supp 0(0) . 

By Theorem El 

<p(t) = <j) f (t) + Mt) , 

l<i<n 

where (f>i(t) denotes the spherical wave 

<tn{t) :=6(t-d l )C(t~d, l )G l . 
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Since the free wave equation has velocity of propagation equal to one, 
we have (here and below t > 0) 

diam(supp <pf(t)) < 2t + diam(S'o) 

and 

diam(supp (f>i(t)) = diam({<ii < t}) = 2t . 
Thus, if Y C S , 

diam(supp <f>(t)) < 2t + diam(S'o) . 

Let us now consider the case Yq := Y \ (So fl Y) ^ 0. 

Suppose that A and B can be chosen diagonal. Let jji G Y . By (TTTfll . 
in order the wave (f>i(t) be present, the point yi must have been reached 
at some earlier stage either by the free wave 0/ or by a spherical wave 
4>j originated from a different point yj. Since all these waves travel 
with the speed v = 1, in conclusion one has that 

diam(supp <f>(t)) <2t + diam(5'o) . 

We have proven sufficiency 

To prove the necessity consider the operator —A A ' B described by 
arbitrary matrices A, B satisfying Q and (J1J. Take an arbitrary point 
yk G Y and initial data satisfying the following conditions: 

4 = &(0) = 0, (21) 

dist^So) <dist(F\{ ?/fe },5o)- (22) 

Then the solution <pf of (JHJ), ([To]) reaches the point before it 
reaches any other point from Y. Therefore, for some sufficiently small 
e > and for any t such that dist(y&, So) < £ < dist(^, So) + e, the 
system of equations (|16|) takes the form 

r ^iCj + SyCj = ^0/(y*) , i = 1, n . (23) 

l<i<" 

The evolution has finite speed of propagation only if all Q, i ^ A; are 
zero, otherwise 

diam(supp <f>(t)) > \y { - y k \ , dist(y fc , 5 ) < t < dist(y k , S ) + e 

which violates (|TT|) since both diam(6 , ) and dist(yk, S ) can be made 
arbitrarily. 

If all C, are zero for % ^ k then ()23|) becomes 

— A ifc C{ + ^ifeCj = A ik (f) f {y k ) , z = 1, n , 

so that the k-th columns of A and B are linearly dependent. Since 
the point y k is chosen arbitrarily from Y, finite speed of propagation 



FINITE SPEED OF PROPAGATION AND LOCAL BOUNDARY CONDITIONS 9 

implies that the linear spaces Vk spanned by the k-th columns of A 
and B are one-dimensional for all k. Then condition (@J) implies that 
any vector v = Ylk=i v ki v k Vk is equal to zero only if all Vf. = 0. 
Therefore the corresponding Lagrangian plane is a tensor product of 
Lagrangian planes for boundary values associated with each singular 
point %)ki i.e. that the boundary conditions are local. 

□ 
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